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Abstract 

The evolution of thermomagnetic perturbations in the resistive state of superconductors is con- 
sidered. A qualitative pattern of the formation and further development of nonlinear stationary 
structures that describe the final stage of thermal and electromagnetic perturbations in a supercon- 
ductor is investigated. The wave propagation velocity and the wave front width in a superconductor 
are estimated. 

Energy dissipation during the motion of vortices leads to ohmic heating of a superconduc- 
tor. As a result, a certain region of the superconductor is heated to a temperature T > T c , 
where T c is the critical temperature. An increase in the temperature in a local region of the 
sample brings about a decrease in the critical current j c and the emergence of vortex electric 
field E at the same region. The effect of superconductivity quenching due to thermal heating 
of the vortex lattice has been experimentally studied for a long time. Early experiments [I] 
revealed that the vortex electric field E is induced under ohmic heating in a superconductor 
through which a direct current with a density j c flows. According to the balance between 
dissipative and nonlinear effects, the transition to the resistive state is accompanied by the 
emergence of various modes of the "switching wave" type, i.e., the regime of a wave motion 
that switches a sample from the superconducting to the normal state. Examples of these 
modes can be provided by thermal waves, namely, a steady-state propagation of a normal 
zone [2] or nonlinear thermomagnetic waves [3] in superconductors. 

In this work, we studied the qualitative pattern of formation and the profile of nonlinear 
dissipative structures, i.e., stationary traveling waves which describe the final stage of the 
evolution of thermal and electromagnetic perturbations in the resistive state of supercon- 
ductors. 

The evolution of thermal and electromagnetic perturbations in a superconductor is de- 
scribed by the nonlinear one-dimensional heat conduction equation 

dT d 2 T 



the Maxwell equation 

d 2 E _ 4vr dj 
dx 2 c 2 dt ' 
and their related equation of the critical state 
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j = j(T,H)+j r {E). (3) 

where v is the heat capacity and k is the thermal conductivity coefficients; j c and j r are the 
densites of the critical and resistive currents, respectively. 

The model under consideration is essentially nonlinear, because the right-hand side of Eq. 
(1) contains a term describing the Joule heat generation in the region of the resistive phase. 
An exact solution to the essentially nonlinear parabolic partial differential equations (l)-(3) 
does not exist. 

Note that the evolution of perturbations of the temperature T(x,t) and the fields E(x,t) 

and H(x,t) is determined by the equation of critical state (3). Because of considerable 

analytical difficulties, we will restrict ourselves to considering the Bean model [4] and assume 

dj c 

that the critical current density is independent of the external magnetic field; i.e., = 0. 
The dependence j c (T) is described by the relationship j c (T) = j — a(T — T ); where T 

is the initial temperature of the superconductor and the quantity a = — - describes 

d-L T=T 

the thermally activated weakening of the Abrikosov vortex pinning by lattice defects. The 
dependence j r {E) in the region of the electric field E > Ej (Ef is the boundary of the linear 
section in the current-voltage characteristic of the superconductor) can be approximated 
by a piecewise-linear function j r {E) ~ <JfE, where o - / is the effective conductivity. The 
dependence j r {E) is essentially nonlinear in the flux creep region E < Ef [5]. Here, we will 
consider a perturbation with a sufficiently high amplitude and use the linear dependence 
ME). 

We seek a solution to the initial system of equations as a function of the new self-simulated 
variable £(x,t), that is 

T = G[£(x, t)],E = E[£(x, t)],H = H[£(x, t)\ , (4) 

which £ = x — v t describes a traveling wave that moves with a constant velocity v along the 
x axis [3]. 

After substitution of relationships (4) into the initial system et of equations and performing 
a trivial differentation, we obtain the following system of equations for the variable £(x,t) 
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+ jE, (5) 



dE _ Attv 
It = 
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E = -H. (7) 

The corresponding thermal and electrodynamic boundary conditions for Eqs. (5)-(7) have 
the form 
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dT 

T(e-+oo) = T , _(^-oo) = 0, (g) 
£(£^+oo) = 0, £(£ -> -00) = £ e ■ 

where _E e is the amplitude of the electric field. It should be noted that inclusion of the 

temperature dependences of the parameters k and v substantially complicates analytical 

calculations of the wave evolution dynamics that is described by the system of Eqs. (5)- 

(7). In most cases, the changes in the local values of these parameters in the sample can 

be considered negligible comparatively to the characteristic scale of temperature variations. 

Hence, we can take these parameters to be constant. Indeed, the investigation revealed 

that the thermal conductivity almost does not affect the character of the stationary wave 

dT 

propagation. This stems from the fact that the thermal flux k — vanishes at stationary 
points of the system at £ — > ±00. However, the temperature dependence of the heat capacity 
should be taken into account. Such a dependence is representeq as v ~ u ( — ) over a wide 



range of temperatures [5]. 

By eliminating variables T(x, t) and H(x, t) with the aid of relationships (5) and (7) and 
employing the boundary conditions (8), we obtain a differential equation for the E wave 
distribution: 



dz 2 n E K c 2 
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Here, we introduced the following dimensionless parameters: 

f _ cH e _ k n_vt K Ana f K _ vL 2 

z — T •> *-* ~ ~a — ~ ' k — —rz , p — —r- , T — — , l K — 



Attjq aL z L c A v k 

where L is the depth of the magnetic field penetration into the sample and t K is the thermal 
time of the problem. According to the qualitative theory [6], the equilibrium states are found 
from the condition 
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E 2 . (10) 



An evident property of system ( 10) is the absence of closed curves that are fully composed 

(dE \ 
E, )■ The proof of this statement can be 

based on the Bendixson criterion [7] . The number of stationary points ( one or three) and 
their type are determined by the parameter 

v 2 

P = 2ttu T —. (11) 

c 2 

The three equilibrium points E — 0, E — E\ and E = E 2 correspond to the condition 
P < Pk = — (Fig.l). There is only one singular point E = at P > P^. The parabola 
and the quartic curve in Eq. (10) are tangent at P = P^, i.e., this condition corresponds to 
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the coincidence E\ — E 2 — E* — The direct solution of Eq. (10) yields the following 
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waves: 



1) E 1 , 2 = E*[l + 2,2(P k -P)^], 
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2)E 1= S^ V -, (12) 



a c 2 



r 2,,2 



E 2 = (2tt) 1 / 2 ^- » , at P fe »P. 



a 2 v 



Analysis of the phase plane shows that the points E = and E = E 2 are stable nodes and 
that E = Ei is a saddle. In addition to the separatrix EiE , the set (10) has the separatrix 
EiE 2 connecting the points Ei and E 2 E\ (Fig. 2). This means that two types of waves with 
amplitudes AE = E 1 and AE = E 2 — E\ can exist in the superconductor. Evidently, wave 
1) has an amplitude of the order E k at P — > P k ; its velocity is determined by equality (10) 
at E — Ei. Equation (10) has two stationary points at P « P k : Eq = is a stable node 
and Ei = 2(3 2 rE K is a saddle. 

The separatrix that connects these two equilibrium states corresponds to a " difference" - 
type solution with amplitude E e , which is related to the wave velocity ve and the wave front 
width Az by the following equations: 
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V E = T 



E P 



2tE, 



(13) 



Wave 2) has a small amplitude at P — > P k 

AE 



E 1 1/2 

. (14) 



E k 



E e 



= 4,4(P fe -P) 1 / 2 «1, (15) 



and its velocity is inversely proportional to the amplitude at P << P k . Such an exotic 
dependence of ve on AE = E e most likely means that the waves of this type are unstable. 
Note that observation of the second-type waves becomes possible in finite-sized samples 
with asymmetric boundary conditions. In conclusion, it should be noted that the above 
investigations prove the possibility of applying the results obtained to high-temperature 
superconductors cooled to liquid-nitrogen temperatures (T = 77 K), providing that the 
values of the physical parameters of the sample are known. 
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